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Abstract. This paper addresses the question of change of decay rate from 
exponential to algebraic for diffusive evolution equations. We show how the 
behaviour of the spectrum of the Dirichlct Laplacian in the two cases yields 
the passage from exponential decay in bounded domains to algebraic decay 
or no decay at all in the case of unbounded domains. It is well known that 
such rates of decay exist: the purpose of this paper is to explain what makes 
the change in decay happen. We also discuss what kind of data is needed to 
obtain various decay rates. 



1. Introduction 

The purpose of this paper it to address the following two questions: 
o What makes solutions to diffusive evolution equations, with underlying linear 
part modelled by the heat equation, dramatically jump from exponential decay, 
when considered on a bounded domain f2 C R™ , to algebraic decay or decay without 
a rate when considered on the whole of W 1 . This will be referred to the decay- change 
phenomenon (DCP). 

o What conditions are required on the data to ensure specific rates of energy 
decay? 

It is well known, that solutions to the heat equation (and solutions to similar 
linear second-order parabolic partial differential equations) defined on a bounded 
Lipschitz domain C M™, n > 1, with initial datum it in L 2 (S1), and subject to 
homogeneous Dirichlet boundary condition, decay exponentially in the L 2 (f2) norm. 
This is an easy consequence of the Poincare inequality: 

V«eHj(0) AxIMI 2 . < ||Vw||l, 
where Ai > is the smallest eigenvalue of the Dirichlet Laplacian on f2, defined 

by 

-A : v e i ► -Av e H _1 (fi). 

It turns out that Ai :— Co/d 2 , where d — diam(£l) is the diameter of f2 and Co is 
a positive constant dependent only on the shape, but not the diameter of f2. This 
is easily seen by performing the change of variable fl C K" i— > Cl = ^(£1 — xq) C R" 
where Xq is the barycenter of f2, and noting that diam(fi) = 1. Thus, as the domain 
grows, lim^oo Ai = 0, and the Poincare inequality is lost in the limit of d — > oo. 
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Unsurprisingly, on R" one can construct solutions to the homogeneous heat 
equation that decay at a very slow algebraic rate, and even ones that decay but 
without any rate. Moreover, we can find, for each time T > and each e£ (0, 1), a 
solution u to the homogeneous heat equation u t — Au = on R" with initial datum 
u e Vp {v : \\v\\2 — (3 < 00} such that 



One can also show the following theorem. 

Theorem 1.1. There exists no function G : (0,oo) 2 1— > G(t,(3) G R+ such that, if 
u is a solution to the homogeneous heat equation u t — Am — on R™ with initial 
datum uq GV/3, then 



Analogous results hold for solutions to many nonlinear equations with a diffu- 
sive term modelled by the Dirichlet Laplacian or a fractional Dirichlet Laplacian, 
including the Navier-Stokes, Navier-Stokes-alpha, the quasi-geostrophic, and the 
magneto- hydrodynamics equations, (see, for example, pQ, [6]). 

To explain the DCP, we shall take a careful look at the behavior of the spectrum 
of the Dirichlet Laplacian near zero on a bounded domain £1 C 1". Although the 
distance to zero of the smallest eigenvalue of the Dirichlet Laplacian on decays 
to zero as the diameter of f2 tends to infinity, for any fixed f2 the distance of the 
smallest eigenvalue to the origin remains positive, so the spectrum, which is entirely 
discrete and 

consists of eigenvalues, remains bounded away from zero. This, as we will try to 
show, is the reason why there is no transition to slower than exponential decay, as 
t — ► 00, of \\u(-, t)\\z for solutions u of the homogeneous heat equation u t — Am = 
in bounded domains Q, subject to homogeneous Dirichlet boundary condition. We 
will clarify this issue by first considering the heat equation in a very simple one- 
dimcnsional situation (viz. on bounded intervals), and then extending the results 
to higher dimensions. For unbounded domains we will derive an extension of the 
Poincare inequality, which will highlight the role played by neighborhoods of zero 
in frequency space. 

This extension of the Poincare inequality will be shown to be optimal in a sense 
that will be explained below. 

The final sections will focus on finding the most general data for solutions to 
the heat equation in unbounded domains so that the corresponding solutions decay 
at a specific rate, then extend these ideas to certain nonlinear evolution equations 
whose solutions satisfy energy inequalities of the type 



\H;t)h<G(t,0) 



and 



lim G(t, /?) 



for all f3 > 0. 



[ \u(x,t)\ 2 dx <-C [ \Vu(x,t)\ 



2 dx. 




2. Notation 



In this section we collect the notation that will be used throughout the paper. 
The Fourier transform of v £ S(M. n ) is defined by 
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u(0 = (2tt)- b f e-' lx ^v(x) dx, 

extended as usual to S'. For a function u : M™ — > C and a multi- index 7 = 
(71, 72, ... , 7 n ) S Nq , -D 7 w denotes differentiation of order |7| = 71 + 72 + • • • + 7« 
with respect to the n (spatial) variables. If v also depends on time t, the symbol Dfv 
is used to denote jth derivative of v with respect to t. If k is a nonnegative integer, 
W fc,p (M") will signify, as is standard, the Sobolev space consisting of functions in 
L P (IR™) whose generalized derivatives up to order k belong to L P (R"), 1 < p < 00. 
When p = 2, W k > 2 (R n ) = H fe (M n ), where the space H S (R") is defined for all s e K 
as the space 

of all feS' such that (1 + |£| 2 ) s/2 /(£) E L 2 (R"). 

3. Preliminaries 

In this section we recall some well-known facts regarding the basis of eigenfunc- 
tions for the Dirichlet Laplacian in bounded Lipschitz domains £1 C R n , and their 
connection to the Poincare inequality. 

The next result is well known (see [3], for example). 

Theorem 3.1. There is a countable orthonormal basis for L 2 (17) which consists of 
eigenf unctions of the Dirichlet Laplacian —A : Hq(Q) H _1 (J7). The eigenfunc- 
tions belong to Hj(f2) and the associated eigenvalues are all positive and bounded 
away from zero. 

To begin, we will suppose that n = 1 and £1 = (—R,R). 

Corollary 3.2. The set of functions, where n = 0, 1, 2, . . . 

„ f [2 , /mrx\ [J f(2n+l)TTx\ 1 

6 = \vr m hr)' Vi cos l 2R )>] 

form an orthonormal basis for L 2 (— R, R). 

Proof. In the light of the previous theorem it suffices to seek A £ C such that the 
two-point boundary-value problem 

d 2 w 

\w + -— = 0, xe(-R,R), w(-R) = w(R) = 0, 
dar 

has a nontrivial solution w G H 1 (— R, R). There is a countable set 

{A„ : n= 1,2,...} U {A„ : n = 0,l,...} 
of A's for which such functions w exist. An easy computation shows that 




sin(2gs) when A = A„ := (2f) 2 , n e {1, 2, . . . }, 

coe^iasg™) when A = A„:= (i 2 ^) 2 , ne{0,l,...}. 



Hence, due to Theorem 13.11 the collection of functions in the statement of the 
Corollary form a basis for L 2 (— R, R). 

□ 
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Let {A„, A„}„ be as defined above. As a consequence of Corollary [321 anv real- 
valued function u € L 2 (—R,R) can be expanded in terms of a countable basis 
consisting of complex exponentials, as follows: 



(3-1) «C»0 = 2i2 



OO OO 

exp ^i-^-ir^ + u(n) exp 



7i— — oo.n^O 



. (2n + l)7r 
1 2i? £ 



Here u(n) = u(—n) and £t(n) = — u(— n) (thus ensuring 
that u(x) € K for all x £ (-R, R)), where u n , n £ Z\ {0}, 
and «„, n e Z, are the corresponding Fourier coefficients. 

The main fact to note is that there is no term corresponding to n = in the first 
sum. 

Let us now take an atomic measure v supported on the set 



(3.2) M = U^} n , ne Z\{0}; 



(2n + 1)tt 
2R 



n € Z 



and let /Lt^ = g^f. Note that rnin^gM |£| = Then, the last expression (|3.1|) 
can be rewritten as 



(3.3) u(a;) = / exp(i£a;) d/i£ = / u(£) exp(i£:r) d/i£. 

J M 

We chose to write u as (13.31) so as to ensure that we have a unified notation, 
regardless of whether we work on — (-R, R) or the whole of K. Hence, by 
Parseval's identity, 



Nli= £ |«W| 2 + 2 

n— — oo,n^0 n— — 00 

Suppose now that t£ G Hq(— 7?, Then, 



■. OO 

d . . v-^ irar v /n7Tx\ 



n— — oo,n^0 



^ i(2n+lV^, , /(2n + l)7rx 
2^ — ^ — exp I 1 



n— — oo 



2R 



R 



Hence the L 2 (— R, R) norm of the derivative is 
(3.4) 



du 
cl.r 



E 



I nir 



i(n)| 2 + £ 



71— — OO 



(2n+ 1)tt 



2i? 



u(n) 



2 n= — oo,ra^0 

On rewriting the last inequality in (|3.4p in integral form, using the atomic mea- 



sure [i£ defined above, we obtain 
(3.5) 



du 
dx 



<w > (^) ii«m. 
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Remark 3.3. The results we obtained above on fi = (-R, R) C R can be easily 
extended to O = (—R,R) n C R n . The eigenfunctions are then products of sines 
and cosines. Hence we can, again, expand a real- valued function u e L 2 (f2) into a 
convergent (in L 2 (fl)) infinite series of 

complex exponentials. Upon doing so, we can express u as the integral with 
respect to an atomic measure ^ supported on the countable set that comprises the 
(discrete) spectrum of the Dirichlet Laplacian on £1 = (— R, R) n . 

(3.6) u(x) = / u(£)exp(i£ • x) d/x^ = / u(£)exp(i£ • x) d[i£. 

JR" Jm 

As in the case of n = 1, the support of the measure ^ excludes £ = 0. 

4. A POINCARE-LIKE INEQUALITY ON R n 

Poincare's inequality is not valid when Q = M. n . We shall nevertheless show 
that a modification of Poincare's inequality holds on R™ by using Fourier transform 
instead of Fourier series 

Theorem 4.1. For each u e H 1 (IR n ) and any A > 0, the following inequality holds: 

\\\7u\\j>A 2 [ \u\ 2 dt- [ (A 2 -|£| 2 )|£| 2 de 

Proof. This follows immediately by PlanchcrcPs identity and by splitting the fre- 
quency domain as 



l n = S U S c where S = {£ : |f | < A}, 



that 



l|V«||l = / lC| 2 |«| 2 de+ / \e\u\ 2 dti 
Js c Js 

■ ( \u\^+ ( |C| 2 |«| 2 de = A 2 / \ufdt- ( (A 2 -|e| 2 )N 2 d£, 

Js c JS Jt" Js 



IS" 

> A 2 



and the conclusion of the theorem follows. 

□ 



Remark 4.2. Suppose that ft is a bounded Lipschitz domain in W 1 with Poincare 
constant Cn, that is, 

. , HVvlla 
Cq := mf 



"6HJ(0) ||V|| 2 

If in the last theorem we choose A = ^,we recover the Poincare inequality from 
the bounded domain = (— R, R), minus an integral in frequency space over an 
interval of 'radius' A = |^| centred at £ = 0: 

l|V^)||^>f^) 2 || W ||^ -f^) 2 / \u\ 2 dt 



\2RJ " " z \2RJ 



« : ICI<27?} 



We now show that Poincare's inequality for the whole space is optimal in the 
following sense: 
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Proposition 4.3. Given any K > and any (3 > 0, if there exist a fi > such 
that 



Nil < j^\\Vu\\l + [x Vu e H 1 ^"), with \\u(x)\\l = (3 



the 



Proof. We show by an example in R 2 that given any \i we can find a function 
u G i? 1 (R 2 ) 3 \\u(x)\\?, — (3 for which the corresponding Ak(u) is smaller then 
[l. Similar examples can be found in any R n only that the computations are more 
tiresome. 

We will work with the family of functions u a = ^|ae _ ' = 1. Straightforward 
computations shows that 



(4.1) Il u a|l2 — P ^ a an( l u a = -^£a 1 e 



v2 I e |2 . 



(4.2) ||Vu Q || 2 =l3a 2 and / \u a (^)\ 2 d£ =(3(1- e [ ^ ] ) 

Jmt\<K} 

Now using (|4.ip and (|4.2p in the Poincarc inequality for the whole space gives 



(4.3) (3< —n + Pil-e^) 

Let f3{l - e [! £ ] ) = B(a) 
We want to show that if 

(4-4) (3<^{3 + n 

then fi > (3(1 - e [J £ ] ) 

When a > K 2 then the first term in the RHS in (|4.3|) is larger that the LHS of 
the inequality so fx — will suffice to have the 

(4.5) p<^ n + fl 

We need to show that 



(4.6) fj, > B(a) = 13(1 - —) 

for all a, in particular if we take a = =|- it follows that it is necessary that /i > 4 
Noting that 

lim B(a) = 

it follows easily that we can find a close to zero so that B(a) < \ 
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□ 

4.1. Poincare for Poisson type equations. We now describe a different way we 
can write a Poincare type inequality for the whole space. We will call this a "Fake 
Poincare inequalit" (FPI) 

Proposition 4.4. Let u e H 1 (R n ) then there exist constants K and a < 1 so that 
(l-a)\\u\\l<K-i\\Vu\\l (FPI) 

J s u d£ 

Proof. Use Poincare for the whole space with a = — jj^p — . Where Sk = {£ : \Q < 
K} and we choose a K so that J s u d£ < \\u\\2- □ 



We now show how (FPI) can be used to obtain a Poincare inequality for solutions 
to appropriate differential equations. 

Example 4.5. Given some constant M, 

let Vt = {.f € £ 2 : \f\ < MM™, for |£| < [3} 
Let to > k. Suppose u satisfies 

(4.7) D k u = f 

(4.8) / G Vv^ n , some to, and j3 

(4.9) lim u = 



then 



ll«lli < 2||v«ni 



Proof. Let S a — : \£\ < a} and M a — M\u>\, where \u\ is the measure of the unit 
sphere. 
Note that 

~ — L 

• U- j^p 

• Is a l"l 2 « = Is a W* d ^ M [3 2 ™-^ 

• Let [3 be such that M o 1 J n ~' lk+n < ^ 

Then if we apply (FPI) the conclusion of the example follows. 

□ 

Remark 4.6. Note that cL 2 n I m , where I m = {/ : I m (f) e ^(K")}, with 7 m 
the Riesz potential of order to for the function /. 

4.2. Poincare and decay. Let u(x,t) G i/ 1 (M™) satisfy 

(4.10) |lK*)lli<-c||Vu||i 

The Poincare and the modified Poincare inequalities can be applied (for bounded 
or unbounded domains respectively) to the RHS of last inequality. 

For bounded domains to express the Fourier series as a Fourier integral, one 
uses a measure which has discrete support. Since zero is not in the support the 
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decay will the modified Poincare inequality for the whole space that will induce the 
passage from exponential to algebraic decay or decay without a rate. 

To get decay on the whole space for solutions u of the inequality ([4.10[) we use 
Fourier splitting. 

This method applied to the whole domain, shows that we can look at frequencies 
near the origin in balls that depend on time. More precisely Fourier splitting is just 
an application of the modified Poincare inequality in the whole space 

where we have chosen A = A(t) appropriately. 

Theorem 4.7. Let u be a solution to 



(4.11) u t =NL(u) + Au 

(4.12) m G L 2 (R n ) 

In addition suppose that the following properties hold 

• J K „ u ■ NL{u) dx = 

• |«o(0l < c\£\ k f° r l£l << !> some k >° 



• NL(u) can be approximated in L 2 by NL(u n ) where u n are sufficiently 
smooth. 

Then 

there exist a constant C = C{u ) so that 

\\u(t)h<C(t + l)-^ 

Proof. Apply Poincare for whole space and then Fourier Splitting to approximating 
sequences of solutions . See [5] for this method. Then pass to the limit. 

□ 

Remark 4.8. We note that the same theorem can be applied in case that diffusion 
is described by (— A) a u with a fractional, or by D m u. The decay rates will have 
to be changes appropriately. 



5. Classification of data 

As seen through Theorem II .![ for arbitrary data just in L 2 there is no specific 
decay rate for solutions of the heat equation. That is for any fixed energy value we 
can find data which leads to a solution whose heat energy decays arbitrarily fast or 
slow. Hence, the energy decay rate is dependent on the actual fom of the data and 
not on initial energy. In this section we characterize the type of data which leads 
to exponential or algebraic decay. 

It is known that the structure of the data near the origin in Fourier space dictates 
the rate of heat energy decay and our theorems rely on this relation. In particular 
we show that exponential decay can occur if and only if the data is zero in some 
neighborhood of the origin in Fourier space, we show that if a Riesz potential of 
the initial data is in L 1 this can determine decay rate, and finally we introduce a 
way to find what type of polynomial "best" describes the data near the origin and 
use this to determine decay rates. The last piece relies of finding a unique decay 
character for any L 2 function by examining the Fourier transform near the origin, 
this can be used to determine the algebraic decay rate exactly. We first will analyze 
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solutions to the Heat equation then extend the results to a more general setting of 
parabolic equations which have a Laplacian linear part. 

The starting idea is that functions with Fourier transform equal to zero near the 
origin decay exponentially, this suggests that a band pass filter will be useful in 
characterizing such functions. 

Lemma 5.1. Given p > 0, let Xp(0 be the cut-off function equal to 1 when < 
pj, j = 1, ...,n and equal to zero elsewhere. A function u G L 2 (M") satisfies 
= for a.e. < pj, j — 1, ...n if and only if u = v — H p * v a.e. for 

some v £ i 2 (R") where H p (x) = rij=i H Pj (xj) = JTj=i slD £ jXj . (Here H is used 
to denote a high pass filter). 

Proof. This proof is quickly checked by noting the Fourier transform of H p * u is 

a 



The following theorem establishes that a solution to the Heat equation decays 
exponentially if and only if its initial data is zero a.e. in some possibly small ball 
centered at the origin in Fourier space. Such functions have a particular form in 
the original space, shown by the previous lemma. This characterizes all data which 
leads to exponential decay of heat energy. 

Theorem 5.2. The solution of the Heat equation with initial data u e L 2 (R n ) 
satisfies the decay bound \\u(t)\\2 < Ce~ ta for some a > and C > if and only 
if the initial data is of the form u = v — H p * v for some v € L 2 (R n ). 

Proof "<^" 

Let B = {£:|&| </>,•} 

Assume u = vq — H p * w , by the previous lemma |u(£)| 2 = for a.e. £ e B. 



M\\l= I e- 2|?|2 *K| 2 ^ + n / e- 2|£|2 *|no| 2 ^ 

< f e- 2 M 2t \u \ 2 dt + e- 2 » 2t ( \u \ 2 dt 

J B J B c 

< I e- 2 ^ 2t \u \ 2 d^ + e- 2 P 2t 
Jb 



\uo\\i 



Where p = min{pj\j = 1, ...n}. By assumption the first integral on the RHS is 
zero. 



By way of contradiction assume there is an a and a C so that ||u(£)ll2 — Ce ta 
and J| ? | <p |m(£)| 2 d£ > for all p > 0. Then, for p = ^, there is a c > so that 
i|f|< |w(£) I 2 ^£ > c - This implies: 
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Ce- ta2 > \\u\\l 

J\i\<p J\t\>p 

>e- 2t P 2 [ u 2 d£ 
J\H\<p 



> ce 

Taking t sufficiently large violates this bound. 

□ 

The next Theorem attempts to characterize types of functions which lead to 
slower then exponential decay. The first piece of this puzzle is the Riesz potential 
of initial data, Ip {uq). It is defined in Fourier variables as 

(//3«o)(£) = -re- 
write, when the limit exists, 

(5.1) Ap(u ) = lim ^ = / {Ipu ){x)dx 

Note that this exists for all Ipuo G L 1 . 

Theorem 5.3. Let u be the solution of the heat equation corresponding to n £ L 2 . 
If Ipuo G L 1 then 

t^\\u(t)\\ 2 <C(A (u o )) 2 

Where Ap(uo) was defined above. If \\u(t)\\2 < C(l + t)% + P for some C and (3 
then 

.. . f uo(0 „ 
hm mf -7-7-3- < 00 

Proof. Both statements in the theorem are consequences of the following equality, 
proved by the change of variables y/t£ = V- 

tf +0 \\u{t)\\t = t% + P f e- 2 ^ 2t u 2 (0d^ 
= f e-^ 2 V 2 "(f^(r,)) 2 d V 

JR™ 

If lim^i^o ifjl^ exists and is finite then the Lcbcsquc dominated convergence 
theorem proves the first statement. If ||u(£)||2 < C(l + t)~% ~$ then 



/ e-^^{Ipw{ji)) 2 dr,<C- 
Jr" ( 



Fatou's lemma then proves the second statement. 

□ 
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The above theorem showns when Ap(uo) exists we can determine the rate of 
decay from (3 and the dimension of space but the condition Ap{ua) exists is stronger 
then i*o € L 2 so a more general structure is needed. Moreover, determining if 
Ap(uo) exists entails determining if Ipuo S L , which might not be always simple. 
On the other hand, if the Fourier transform of initial data is a polynomial of the 
form in some possibly small neighborhood of the origin then the solution will 
decay as C(l + t)~ q ~z . This can be checked by calculating (or estimating) the 
integral J B r p \ e~ 2 '^' *|£| 2l2 <if\ Unfortunately only a small amount of initial data can 
be described in this way and we are pushed to find what order of polynomial best 
describes a general function uq 6 L 2 near the origin in Fourier space and using this 
information to determine the decay rate. 

Definition 5.4. The "decay indicator" P q (uo) is defined as follows. Let B(p) be 
the ball of radius p centered at the origin, for q € (— oo): 

P q (u ) = limp- 2 '-" / \uo(®\ 2 <% 
Jb( p ) 

When m £ L 2 (R") the integral j B i p \ \uo{£)\ 2 d£ considered as a function of p 
is continuous, monotone decreasing, and bounded from below as p becomes small, 
this is enough to ensure that P q is defined for all q and u n e L 2 . It is trivially zero 
for all q < — § when uo G L 2 so we consider only q S (— §,oo). P q (u a ) compares 
uq to the polynomial |£| 2? near the origin and takes values in the (non- negative) 
extended real numbers, we are interested in three outcomes: P q {uo) — 0, oo, c with 
c 7^ 0. Depending on the outcome we think of uq as a polynomial with degree, 
respectively, greater, less, or equal to q near the origin. It is easy to check that 
P q (m = \L0\(n + 2q) = p q . 

Recalling (|5.ip . when A q (uo) exists it bounds P q (uo), in this sense P q is weaker 
then A q . Indeed, 

P q (u )= limp- 2 ?-" f \u \ 2 d£ 
Jb(p) 

< lini J sup p- 2 ?- / |C| 2 "^ 

\b(p) iei 2? / Jb(p) 

Definition 5.5. For a given uq S L 2 (R"), we call the unique value q* given by 
Lemma WM below the "decay character." 

Lemma 5.6. For each uq G L 2 there is at most one value of q € (— §,oo) so that 
< P q (uo) < oo. If such a number exists we denote it q* . If no such number exits 
we take q* = — ^ in the case where P q {uo) = for all q S (— §, oo) and q* — oo if 
Pq(uo) = oo for all q £ (— f,oo). 

Proof. Let a = sup{<7 : P 9 (uo) = 0} and b = inf{q : P q (uo) = oo}. If q is such that 
P q (uo) — oo and r > q then P r (uo) — oo, this is observed by taking the limit of the 
following inequality which is valid for all p < 1: 
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(5.2) p- 2r ~ n f \MO\ 2 d( = P 2q - 2r p- 2q - n f \MO\ 2 dt 

Jb(p) Jb(p) 

> P - 2q - n l \M0\ 2 <% 

From this we conclude P q (uo) = oo for all q > b and the lemma is true if b = — ~, 
a similar statement can be made concerning q < a and shows the lemma is true if 
a = oo. 

It is also clear that a < b and q* will be well defined for all uo E L 2 when 
the proof is finished. To accomplish this we wish to show a — b, this is proved 
by contradiction. Assume, contrary to the statement, there exists q £ (a, b), then 
Pq(uo) = c for some < c < oo. If e > 0, similar to the above inequality: 



Pq+e{uo) = 




This shows P 9 + c (ito) = oo, since e was chosen arbitrarily we conclude q = b. 

□ 

The decay character is calculated from the behavior of an L 2 function near the 
origin in Fourier space, we now prove a theorem relating decay rates and P q (uo). 
A consequence of this theorem (Theorem 15.81) will summarize the relation between 
the decay character and the exact decay rate which is our goal. 

Theorem 5.7. Let u be the solution to the heat equation corresponding to u € 
L (W 1 ) and q £ (— ^,00). // P q (u ) > there exists a constant G\ > which 
depends only on HmoII! an d the dimension of space so that 

Ci(i + <II«(*)II2 

If P q (u ) < 00 there are constants C2,C% > 0, again depending only on || u 1 1 1 
the dimension of space so that 

\\u(t)\\l<C 2 (C 3 + 

Proof. We consider first the lower bound assuming P q (uo) > 0. Relying on the 
existence of the limit and the fact that it is bounded away from zero we may take 
po > sufficiently small to ensure, for all p < pa : 

Cl <p- 2q - n [ \u \ 2 d(; 
Jb(p) 

Let < p(t) < po, we will soon chose it exactly. 

IK*)II5= I e~ 2 ^ 2t \uo\ 2 dZ+ [ e- 2 ^ 2t \u \ 2 dZ 

Jb(p) JBC{p) 
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Choosing p(t) = po(l + t)~^ proves the lower bound. 

To prove the upper bound assume P q (uo) < oo and take po > sufficiently small 
so that for all p < po'. 

p- 2q - n I \u»?dt<c 2 

Jb(p) 

The constant c 2 is known to exist since P q (u ) < oo. We now use the Fourier 
Splitting Method ([5]), starting with the well known energy inequality for the heat 
equation with < p(t) < p . 

~IK*)lla<-||V«(*)||l 

<-p 2 [ \u(t)\ 2 d^ 

This implies 

±\\u(t)\\ 2 2 + 2p 2 \\u(t)\\ 2 2 <2p 2 [ \u(t)\ 2 d( 
< 2p 2+2q+n C 2 

Take m > q + ^ and choose p(t) — 2 (Ci+t) 1 / 2 wnere C3 > is large enough to 
guarantee p(0) < po. Solve the differential inequality with the integrating factoi- 
ds + t) m to find 

\\u(t)\\l < C(C 3 + + (C 3 + t)- m \\u \\ 2 2 

This is the upper bound. 

□ 

Theorem 5.8. Let uq G L 2 (M. n ), u{t) the corresponding solution to the heat equa- 
tion, and q* the decay character. If < q* < 00 then there are constants 
Ci, C 2 , C 3 > so that 

(5.3) Ci(l + t)- q "~9 < \\ u (t)\\ 2 < C 2 (C 3 + t)- q "~9 

Moreover, if q* = — ^ or q* = 00 then ||u(i)||2 decays, respectively, slower or 
faster then {l + t)- q ~%. 

Proof. If — t- < q* < 00 then < P q * (uq) < 00 and (|5.3|) follows from the previous 
theorem. If g* = — § we have, for all q € ( — §, 00) , by the previous theorem we are 
guaranteed the existence of a constant C\ (q) > such that 

Ci(i + <|K*)||2 

Letting q take all values in (—^,00) proves the slow decay statement. The 
statement concerning q* = 00 is argued similarly. 

□ 
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6. Decay for solutions to some nonlinear evolution equations 

Our next goal is to apply this idea to a more general class of PDEs with a 
non-linear term. The idea is that if the non-linear term decays faster then the 
linear term, the decay of solutions will be given through the decay character. If the 
non-linear term decays slower then it will determine an upper bound on the decay 
rate, a lower bound will require more knowledge of the specific non-linear structure. 
First we will assume a bound on the non-linear term in Fourier space and derive 
a bound on energy decay. After this we will use results from [5] with the decay 
character to describe energy decay of solutions for the Navier-Stokes equation. 

Theorem 6.1. Let u be a solution of the PDE \^.ll^ and assume the non-linearity 
satisfies the following conditions: 

1. We are justified in writing the solution as 

u = e-l«l 2t iT + f e-l«l 2 (*- s )iVX(u)(£, s) ds 
Jo 

2. NL^)(C,s)<C\^\ k 

3. J R „ u ■ NL(u) dx = 

Let q* be the decay character associated with uq, then for any q < q* there exist 
constants C\ , C-x > so the solution of the PDE satisfies the energy decay estimate 

\\u(t)\\ 2 2 < d{C 2 + i)" 9 "? + d(C 2 + t)- k ~% 

Remark 6.2. The assumptions in this theorem are reasonable for equations such 
as the Navier-Stokes equation, the Navier-Stokes-a equation and the Magneto- 
Hydrodynamic equation among others (see, for example, [T], [5], and [7]). The 
decay rate in this theorem is determined by q* or k, whichever is smaller. 

Proof. Assumptions 1 and 2 imply 

\u\ < \e~M 2t u \+C [ t e-M'<t-')\t\ k ds 
Jo 

Completing the integral on the RHS then squaring yields 

\u\ 2 <e- 2 ^ H \uof + C\e k - 2 
Consider now any q < q* and take po > sufficiently small so that for all p < po 

p- 2q - n I \u Q \ 2 di<C-, 
Jb(p) 

Here C3 is some constant known to exist since q < q* . Assumption 3 allows an 
energy inequality from which to use the Fourier Splitting Method, proceeding now 
as in the proof of Theorem 15.71 

±\\u(t)\\i + 2p 2 \\u(t)\\i<2p2 [ \u(t)\ 2 dH 

at J B(p) 

<p 2 f \u \ 2 di + Cp 2 f |£| fe - 2 d£ 
Jb(p) Jb(p) 

< 2p 2+2q+n C 3 + Cp 2k+n 
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Take m > max{q + f , ^ + §} and choose p(t) = 2 (cJ+t) 1 / 2 wnere C2 > is large 
enough to guarantee p(0) < po- To finish the Fourier Splitting Method multiply 
the above equation by (C 2 + t) m and solve the differential inequality to find 

< C(C 2 + *)-«-* + C(C 2 + t)- k -% + (C 2 + t)- m \\uo\\ 2 2 
This finishes the proof. 

□ 

In the specific case of the Navier-Stokes equation there has been significant inves- 
tigation into the rate at which a solution approaches a solution of the heat equation 
with the same initial data, see [2] , [4] , [8] . We will now demonstrate how the results 
of this section fit with the main theorem in [8]. 

Theorem 6.3. (Wiegner) 

Let n > 2 be a weak solution of the Navier-Stokes equation on K" , 2 < n < A, 



(6.1) d t u- Au + u- Wu + Wp = Q 

V • u = u(0) = u 
which satisfies the energy inequality 



2 



\\u(t)M + 2 J ||Vu(r)||^dr<||«( fl )||^ 

If \\e At u4l < C(l + t)- a , then \\u(t) - e At u a \\l < K{t){\ + t)- d with d = 
5 + 1 — 2 max{l — a, 0} and 

fe(t) for ; a = 0, with e(t) \ for t — > 00 

C7n 2 (i + e) /or a = 1 

C /or a ^0,1 

Remark 6.4. The actual theorem proved by Wiegner is more general as it includes 
a forcing term, forcing terms are outside the scope of this paper. 

Proof. See [8] 

□ 



Theorem 6.5. Let u be as in Theorem ] 6. 3\ and q* be the decay character associated 
with uq. 

1. If — % < q* < 1 — % then there are constants C\, C2, C3 > so that 

d(l + t)- q '-* < \\u(t)\\ 2 2 < C 2 (C 3 + t)-«"-* 

2. Ifq* > 1- § then \\u(t)\\l < C(l + <)-f-' 3 w/iere /3 = min(g*,l). 

3. If q* = —5 and n = 3,4 then ||it(t)||2 decays slower then any polynomial. 



Proof. This is a combination of the the above theorem and Theorem 15.81 In case 
1, Corollarv l5.8l gives 



d(l + t)---* < \\e^u \\ z 2 < C 2 (C 3 + ty 

n j 3n 

2 ' a ~~ 2 

\\u(t)-e At uo\\ 2 2 <C(l+t)- S ?+ 1 - 2 i* 



while Theorem [O] allows (a = q* + § , d = \ - 1 - 2q*): 
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Notice that the second decays faster. Combining these with the triangle inequal- 
ity proves 1. 

Case 2 is similar, this time Theorem 15.81 gives 

ci(i + t)-*-«* < \\e At u \\ 2 2 <C 2 (C 3 + ty^ 
while Theorem IQl allows (a = f + g*,d=| + l): 

\\u(t) - e At u Q \\ 2 2 < h a {t){l+t)-^- 1 
An application of the triangle inequality shows 

\\u(t)\\ 2 2 < h a (t)(l + t)-^' 1 + C 2 (C 3 + t)-t-«* 
Case 3 also follows from Theorems 15.81 and 

16.31 and the triangle inequality. With a = in Theorem 16.31 we have \\u(t) — 
e A *tto||| < e(t)(l+t)~% +1 while ||e At uo||2 decays slower then any polynomial (recall 
e(i) \ 0). When n — 3,4 the first decays faster, thus, for large t, 

2\\e At u Q \\ 2 < \\\e At u \\ 2 2 - \\u(t) - e At u g\ < \\u(t)\\ 2 

□ 
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